Resilience to damage of graphs with degree correlations 
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The existence or not of a percolation threshold on power law correlated graphs is a fundamental 
question for which a general criterion is lacking. In this work we investigate the problems of site and 
bond percolation on graphs with degree correlations and their connection with spreading phenomena. 
We obtain some general expressions that allow the computation of the transition thresholds or their 
bounds. Using these results we study the effects of assortative and disassortative correlations on 
the resilience to damage of networks. 
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The graphs representing many real networks are char- 
acterized by power law degree distributions 0. The ori- 
gin of these power laws can be traced back to the growing 
nature of real networks and to some effective preferential 
attachment mechanism. This later mechanism implies 
that when new vertices are added to the graph they are 
more likely linked to existing vertices with large degrees 
0, • Recently, there has been a great interest in the 
study of processes running on top of these graphs due to 
their social, technological and scientific relevance. Per- 
colation processes spreading phenomena [1,|EI3j 
the Ising model || and searching techniques 0] are some 
examples for which analytical solutions have been found 
in random graphs with the only constraint given by the 
degree distribution. One of the fundamental results is 
that the threshold characterizing the percolation tran- 
sition or an epidemic outbreak, depends on the ratio 
(d 2 ) lid) of the first two moments of the degree distribu- 
tion 0,0,0, Hence, if (<i 2 ) diverges when increasing 
the graph size, there is no transition in the thermody- 
namic limit. 

The topology of real networks is also characterized by 
degree correlations 0, and, therefore, the extension 
of previous results for uncorrelated graphs is of utmost 
importance. Moreover, it has been shown that growing 
network models with 0] and without 0] preferential 
attachment lead to non-trivial degree correlations. The 
study of models on graphs with degree correlations is 
quite recent [ill H3 . fljl H(| . Some expressions for the 
size of the giant component and related quantities have 
been obtained in Ref. |llj whereas an equation for the 
epidemic threshold has been provided in [J^j. General 
statistical mechanics approaches for models on correlated 
graphs has also been developed in Refs. 0,0- How- 
ever, in contrast to the case of uncorrelated graphs, a 
general criterion for the existence or not of a transition 
threshold has not been proposed yet. A first step in this 
direction has been taken in Ref. yjj for a disease spread- 
ing model. 

In this paper we study the resilience to damage (vertex 



or edge removal) of random graphs with arbitrary degree 
distributions and correlations by addressing the problem 
of dilute (site or bond) percolation on these graphs. We 
report a general equation for the threshold and bound it. 
Besides, we analyze the effect of correlations considering 
some examples of uncorrelated, assortative and disassor- 
tative correlated graphs or their mixture. We conclude 
that assortative correlations can make graphs quite ro- 
bust, even with a finite (d 2 )- On the contrary, disassor- 
tative correlations can make graphs fragile, even with a 
divergent (d 2 ). 

Let us start by considering the set of undirected graphs 
with N vertices and arbitrary degree distribution pd- Fol- 
lowing one end of a randomly chosen edge, we will find a 
vertex of degree d with probability q d = dpdj (d). We fur- 
ther assume correlations between adjacent vertices: The 
conditional probability p(d'\d) that a vertex of degree d' 
is reached following any edge coming from a vertex of de- 
gree d explicitly depends on both d and d' . Consistency 
with the degree distribution requires Yld' P(d'\d) = 1. 
Besides, the joint probability p(d'\d)qd that the two ver- 
tices at either end of a randomly chosen edge have degrees 
d and d! must be symmetric. For uncorrelated networks 
p(d'\d) = qd> that is independent of d. 

The problem of percolation on graphs with degree cor- 
relations has been recently studied 01 usm g the gen- 
erating function formalism. Alternatively, one can use 
a more general statistical mechanics approach [l6j| . In 
this case the size of the giant component, the fraction of 
nodes in the largest cluster, is given by 



S=l-J2Pd(u d ) d , 



u d = }^p(d'\d)(u d > 

d' 



(1) 



(2) 



where Ud is the average probability that an edge con- 
nected to a vertex of degree d leads to another vertex 
that does not belong to the giant component [ll| . 
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Let us generalize this result to the site percolation 
problem. In this case a fraction / of the nodes is re- 
moved from the graph and the new giant component is 
computed. Since the node removal is independent of the 
node degree this is equivalent to replace the original de- 
gree distribution and correlations by: (i) the probability 
that a node selected at random has degree d and it has 
not been removed, and (ii) the probability that if we se- 
lect a node at random and follow one of its edges we end 
in a node with degree d! that has not been removed, i.e. 

Pd^(l-f)Pd, p(d'\d)^(l-f)p(d'\d), (3) 

Substitutin; 



on a simple topological measure, whether or not a given 
graph is robust under vertex removal. 

In the bond percolation problem a fraction / of the 
edges is removed from the graph and the new giant com- 
ponent is computed. Since the edge removal is made 
at random this is equivalent to keep the original degree 
distribution and replace the degree correlations by the 
probability that if we select a node at random and follow 
one of its edges we end in a node with degree d! that has 
not been removed, i.e. 

Pd^Pd, p(d'\d)^(l-f)p(d'\d). (11) 



g Eqs. © in Eqs. JTJ and 101 we get Substitution of Eq. (JTTJ in Eqs. ||TJ and © yields 

S=l~ /-(I-/) £^W> (4) S = l-J>M d , (12) 



u d = / + (1 - f^Pid'WM*'- 1 , (5) 

d' 

where the term — / (/) in Eq. Q ( Eq. J5jl) gives the 
probability of hitting a removed node. One solution to 
these equations is u d = 1 yielding S = 0. This solution 
is valid whenever the equation for the Ud is stable un- 
der successive approximations. That is, if we start with 
Ud(n) = 1 — Pd(n) and compute the successive approxi- 
mation pd(n + 1) then we should obtain that pain) — * 
in the limit t — > oo. For pd(n) <C 1 the last equation is 
approximated by the linear map 

p d {n + 1) = Ldd'pd' (n), (6) 

d> 

with 

Ldd' = (l-f)C dd >, C dd , = (d'-l)p(d'\d). (7) 

The stability of the solution Ud — 1 is then related to 
the largest eigenvalue of Ldd 1 ■ If it is smaller (larger) than 
1 the solution is stable (unstable). Since Ldd' is linear in 
/ the stability condition can be written as 

/ > f c , (1 - f c )A max = 1, (8) 

where A max is the largest eigenvalue of Cdd' provided 
that h-max > 1. If A max < 1 the graph does not have a 
giant component even for / = 0. Moreover, since Cdd' 
is a positive matrix then A max has the lower and upper 
bounds minded' Cdd< and max d ^d' Cdd', yielding 

min <i ( d )nn ( d ) < 1 + A max < max d (d) „„ (d) . (9) 

where 

(d) nn (d) = J2p( d w> ( 10 ) 

d' 

is the average degree among the neighbors of a node with 
degree d 10]. Eq. © can be used to determine, based 



u d = f + (1 - f)J2p(d'\d)M d '^- (13) 

d' 

Note that the only difference between the site and bond 
percolation problems (see Eqs. I4I5|I is the equation for 
the giant component while that for u d is identical. Hence, 
Eqs. © and arc also valid for the bond percolation 
problem. 

In what follows we consider some particular graphs in 
order to analyze the effects of correlations. Depending 
on the monotony of (d) nn (d) the degree correlations can 
be classified in: uncorrelated if it is independent of d, 
assortative or positive if it increases with increasing d and 
disassortative or negative if it decreases with decrea sing 
d. A similar definition has been introduced in Ref. [llj 
using a correlation coefficient. 

For random graphs with no constraint other than the 
one imposed by the degree distribution we have p(d' , d) — 
q d i . In this case the lower and upper bounds in Eq. @ 
are equal giving for the largest eigenvalue 

K% = V ~ 2 ' (M) 

Alternatively one can compute A directly from the eigen- 
value problem of Cdd'- Then from Eq. JHJ we obtain 
1 - / c = l/((d 2 ) / (d) - 2) 4]. Hence, if the second 
moment (d 2 ) diverges the threshold equals 1, i.e. the 
network is robust under random vertex or edge removal. 
Furthermore, consider the case in which the degree cor- 
relations can be decomposed into two components 

p(d'\d) = aq d , + (1 - a)5p(d'\d) (15) 

with < a < 1 and 5p(d'\d) > for all (d, d'). Vary- 
ing the parameter a one interpolates between the un- 
correlated graphs (a — 1) and a graph with arbitrary 
degree correlations given by 5p(d'\d). In this case from 
Eq. © we obtain A max > aA™™ and, therefore, if the 



3 




FIG. 1: Size of the giant component for a graph with pd = 
cd~ 3 ' 5 (2 < d < dmax, dmax = 100) and degree correlations 
p(d'\d) — aSdd' + (1 — a )ld' with a, as computed from Eq. 
@. The dashed line marks the percolation threshold ob- 
tained using perturbation theory (Eq. 11811 ). The inset shows 
the largest eigenvalue relative to dmax as a function of a. The 
points were computed numerically and the line is the pertur- 
bation theory dependency A max / 'dmax = a. 



network is robust for the uncorrelated case it will also 
be robust for any a > 0. This immediately implies that 
any graph with a divergent second moment and a finite 
amount of random mixing of edges does not have a per- 
colation threshold. 

Assortative correlations allow us to show that the di- 
vergence of the second moment is not a necessary condi- 
tion for the absence of the threshold. Let us consider a 
network with degree correlations 



p(d'\d) = a5 dd , + (1 - a)Sp(d'\d), 



(16) 



with < a < 1 and Sp(d'\d) > for all (d, d'). a = 1 
corresponds to a fully assortative graph made up of sub- 
graphs with fixed degree. In this case Cdd< = d'Sdd> (see 
Eq. [7J is already diagonal. The largest eigenvalue is 
A max = dmax, where d ma x is the largest degree. If d max 
diverges for N — > oo then f c = l. For the more general 
case < a < 1 we compute the largest eigenvalue using 
perturbation theory |l8| around a — 1, obtaining 



A* 



,{a) = ad max + (1 - a)C dmax d„ 



(17) 



This result is valid whenever (1 — cx)Cd max d masc *C ad max . 
In general Cd ma!C d ma!C decreases with increasing dmax, re- 
sulting 



^max(a) « adr, 



(18) 



for d max ^ I /a. Hence, for any a > and any un- 
bounded degree distribution we have f c = 1, i.e. there is 
no percolation threshold. In Fig. \l\we show the validity 
of the perturbation theory for a particular perturbation 
8p(d'\d). Thus, as in the fully assortative case, if a > 



and dmax diverges f c = 1. Therefore, we can conclude 
that the divergence of the second moment is not a nec- 
essary condition. 

Let us now analyze if the divergence of the second mo- 
ment is a sufficient condition for f c = 1, using an example 
of a disassortative graph. Consider a vertex with degree 
d > 1 and an edge incident to it. Then with probability 
gd a vertex at the other end is chosen at random among 
all vertices with degree d' > 1, otherwise it is connected 
to a vertex with d! = 1 chosen at random, i.e. 



p(d'\d) 



(1 - g d >)d'p d > 



+ (l-g d )5 d >,iQ(d-l 
gd'd'pd 



e(d' - i)s d , 



9d 



E s 9ssp. 



e(d'-l)9(d-l). (19) 



where Q(x) is the unitary step function (Q(x) — for 
x < and 9(x) = 1 for x > 0). Moreover, the fraction 
of nodes with degree 1 is obtained self-consistently from 
the condition p\ = J2 d >i(^ ~ 9d)dpd- The average degree 
of the neighbors of a node with d > 1 is given by 



(d) nn 



9d 



Y.d>>i9d>d l2 p d > 
J2s 9ssp s 



(20) 



and, therefore, these graphs are disassortative for any 
monotonic decreasing function gd- To analyze the perco- 
lation properties of this graph we computed exactly the 
largest eigenvalue of Cdd' = (d' — l)p(d'\d), resulting 



J2d9 2 d(d-l)dp d 
9s sp s 



(21) 



Hence, the conditions for the existence of a giant compo- 
nent (Amax > 1) or resilience to damage (A max = oo) are 
modulated by gd and, therefore, the disassortative corre- 
lations given by g d have a great impact on the percola- 
tion properties. For instance, let us consider gd = d~ a 
and a power law degree distribution p d = cd f with 
7 < 3 ((cP) = oo). From Eq. it follows that 

(d) n — 1 ~ dr a so that when increasing a the graph 
gets more and more disassortative. Moreover, A ma x di- 
verges for a < a c — (3 — "f)/2 and it is finite otherwise. 
Thus, for small values of a the graph is robust but for 
a > a c it becomes fragile. It is worth noticing that the 
value of a above which the giant component disappears 
(Amax < 1) is larger than a c . Besides, for large degrees, 
the degree distribution of the vertices in the giant compo- 
nent is still a power law, but it decays slower than that of 
the whole graph. Thus, disassortative correlations com- 
petes against the formation of the giant component and, 
the divergence of (<i 2 ) is not a sufficient condition to get 
a robust graph with f c = 1 . 

The connection between percolation theory and mod- 
els of epidemic spreading is well known Two general 
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classes of epidemiological models can be related to perco- 
lation problems, the Susceptible-Infected-Removed (SIR) 
and the Susceptible-Infected-Susceptible (SIS) classes. 
The SIR model assumes that individuals can exist in 
three classes and that once they get infected they can 
not catch the infection again. This model can be mapped 
into a bond percolation problem taking / as the probabil- 
ity that the disease will be transmitted from one node to 
another and the size of the giant component as the size of 
the outbreak. Hence, all the conclusions drawn above for 
the bond percolation problem can be translated to the 
language of epidemic spreading for the SIR model on top 
of graphs with degree correlations, extending in this way 
previous studies in Refs. 0,0] for uncorrelated graphs. 

On the other hand, the SIS model allows individuals 
to move through the cycle of infection so that the preva- 
lence (number of infected individuals) attains a station- 
ary value. The SIS model on top of graphs with degree 
correlations has been recently analyzed in Refs. |l4lll7j. 
They obtained the epidemic threshold (the value of A 
above which the solution with zero prevalence is unsta- 
ble) A = 1/A' max , where A' max is the largest eigenvalue of 
the matrix C' dd , = dp(d'\d). This approach is quite sim- 
ilar to the one presented here for site percolation with 
the remark that C' dd , is different (see Eq. (JJJ)). In fact, 
if yd is an eigenvector of C' dd , — dp(d'\d) correspond- 
ing to the eigenvalue A' then yd/d is an eigenvector of 
C dd , — d'p(d'\d) corresponding to the same eigenvalue. 
This last matrix is that of Eq. J7J), but replacing d' 
by d' — 1. However, this subtle difference makes the 
SIS and dilute percolation different. We have computed 
the largest eigenvalue of C dd , for the disassortative graph 
considered above(Eq. (fP3)l h Taking the limit (gP) ^ 1 
one gets 



AL 



Ed(l - 9d)d?Pd 
E^ 1 -9s)sp s ' 



(22) 



where gd is again a decreasing function of d. In this 
case, independent of the form of gd, the divergence of 
the second moment of the degree distribution implies the 
divergence of h! max . Moreover, the same conclusion is 
obtained if gd is an increasing function of d. The con- 
ditions for the existence of a finite prevalence in the SIS 
model have been recently addressed in ^t|j where the 
divergence of the second moment has been shown to be 
a sufficient condition for the absence of the phase tran- 
sition in the SIS model. Nevertheless, we have shown 
that this conclusion does not hold for dilute percolation. 
This essential difference is rooted in the existence of an 
additional dimension in the SIS model, given by the time 
evolution of the density of infected sites. 

In summary, we have studied the percolation problem 
on top of random networks with arbitrary degree distri- 
bution and correlations, making its generalization to site 
and bond percolation. The connection with spreading 
phenomena was also analyzed. We provide some general 



expressions to obtain or bound the transition threshold. 
Using these results we have shown that the existence of 
a finite amount of random mixing of the connections be- 
tween vertices is sufficient to make the graph robust un- 
der vertex or edge removal provided (d 2 ) — > oo. Assor- 
tative correlations makes the situation even better, they 
can lead to a graph robust to random damage even with 
a finite second moment of the degree distribution. On 
the contrary, disassortative correlations compete again 
the formation of the giant component and can make a 
graph fragile even with a divergent second moment. 
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tially supported by the European commission FET Open 
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